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ABSTRACT 


This thesis addresses two problems in aligning the 
recruiting structure for the Navy Recruiting Command. The 
first problem involves two decisions affecting recruiting 
Stations within a single recruiting district: which stations 
should remain open and how many recruiters should be assigned 
to each open station? The second problem is to decide how 
many recruiters and stations each district should have. 

The first problem is formulated as a nonlinear mixed 
integer programming problem. To obtain a solution with 
readily available software, the problem is decomposed into 
four subproblems that are solved sequentially. This 
decomposition approach is empirically shown to yield near 
optimal solutions for problems of varied sizes. The second 
problem is formulated as a nonlinear resource allocation 
problem in which the objective function is not expressible in 
closed form. To efficiently solve this problem, the function 
is approximated in a piecewise linear fashion using the 
results from the first problem. To illustrate the 
applications of these optimization models, solutions were 
obtained for Navy Recruiting District Boston and Navy 
Recruiting Area 1, which consists of Albany, Boston, Buffalo, 
New York, Harrisburg, Philadelphia, Pittsburgh and New Jersey 


districts. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases of 
interest. While effort has been made, within the time 
available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 


additional verification is at the risk of the user. 
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I. INTRODUCTION 


A. BACKGROUND 

Although the U.S. Navy is becoming smaller, the Navy 
Recruiting Command must continue to recruit in order to 
replace sailors at the lowest grade who are promoted or lost 
to attrition. Military downsizing means fewer recruits are 
needed, but these recruits are more difficult to obtain. This 
is due to several factors. First is the decrease in the 
recruiting budget and the rise in operating costs. Second is 
the emphasis on obtaining a greater proportion of quality 
recruits who are capable of manning ever more sophisticated 
weapons, sensory and communications systems. Recruiting these 
quality individuals is costly because they are in demand from 
other branches of the military, private industry and academia. 
Finally, the last factor is the expected competition from the 
National Service Corps that President Clinton promoted during 
his campaign. Many individuals in the prime target population 
for recruiting would likely find this program attractive, for 
it would be less restrictive and involve less risk than 
military service. In addition, it would keep people closer to 
families and friends, and would offer attractive educational 


benefits. 





To assist Navy Recruiting Command in accomplishing this 
increasingly difficult task, this thesis addresses the problem 
of aligning recruiting stations and districts. Proper 
alignment of these organizational units would allow the 


Command to recruit more effectively. 


B. CURRENT ORGANIZATION 

The Navy Recruiting Command (NRC), located in Arlington, 
Virginia, is the headquarters for a nationwide network of 
recruiters. The command is organized into five recruiting 
areas with each area subdivided into districts. There is a 
total of 41 districts. Each of which is organized into zones 
which are composed of stations. The stations serve as the 
base from which the recruiters operate. In 1991, there were 
1,283 Navy Recruiting Stations in the Continental United 
States. Recruiters actively pursue prospective recruits in 
the territories of their assigned station. In general, the 
territory of a recruiting station consists of a collection of 


adjacent zip code areas. 


C. PROBLEM STATEMENTS 

The problem currently faced by NRC is how to maintain an 
effective and efficient recruiting structure with a smaller 
budget. In this thesis, two problems are considered. One 


concerns the recruiting stations within the territory of a 





single district and the other concerns the allocation of 


recruiting resources to several districts. 

The first problem is to determine, for a given district, 
which stations from a list of candidates to open and the 
number of recruiters to assign to each open station. This 
problem is formulated as a nonlinear integer programming 
problem and it is called the Location-Allocation Problem. As 
is well known in operations research, these problems are 
difficult to solve. However, this thesis obtains near optimal 
solutions to the problem via a decomposition technique. 

The second problem involves the allocation of recruiting 
resources to several districts. For example, if the current 
operating budget only supports 500 stations and 1000 
recruiters, then how many stations and recruiters should be 
assigned to each district? This problem is formulated as a 
nonlinear resource allocation problem and it is called the 
Multidistrict Allocation Problem. The objective function of 
this problem is not expressible in closed form. To solve it, 
this thesis provides a methodology to obtain solutions through 


approximation using piecewise linear functions. 


D. THESIS OUTLINE 

The two optimization models mentioned above use results 
from forecasting models developed by analysts at NRC. rer 
completeness, these models are described in Chapter II. 


Chapter III describes and formulates the location-allocation 





problem. Chapter IV discusses the technique to decompose the 


problem in Chapter III. Chapter V provides a bounding 
technique which validates the claim that the decomposition 
technique nvscvides near optimal solutions. Chapter VI 
describes the multidistrict allocation problem and presents a 
sciution approach. Finally, Chapter VII summarizes the thesis 


and suggests areas for future research. 





II. FORECASTING ENLISTMENT CONTRACTS 


An extremely important input required by optimization 
models for locating recruiting stations and allocating 
recruiters is the ability to forecast the distribution of the 
target market across the continental United States (CONUS). 
Despite the increasing role of women in the military, the 
prime target market for military service continues to be males 
between 17 and 21 years of age with no prior military service. 
From this market, the Navy is primarily interested in 
enlisting quality recruits. The first section in this chapter 
presents NRC’s definition of ‘quality’. The next section 
describes four statistical models developed by NRC to forecast 
the distribution of quality contracts across CONUS. Finally, 
the last section discusses mathematical implications of these 
statistical models when embedded in the optimization models in 
the following chapter. These optimization models work equally 
well with different statistical models. The NRC models are 
selected mainly because of their accessibility and familiarity 


to NRC analysts. 


A. QUALITY CONTRACTS 
The Navy groups prospective recruits into major categories 
or cells based on two factors: educational attainment and 


mental aptitude. For the educational factor, there are only 





two classifications. The first consists of prospective 


recruits who already possess a high school diploma or are 
currently high school seniors expecting to graduate. The 
second consists of those who are currently without a diploma 
and do not expect to obtain one in the future. The mental 
aptitude classification is determined by the score from the 
Armed Forces Vocational Aptitude Battery (ASVAB) or Armed 
Forces Qualification Test (AFQT). Each prospective recruit 
must take the ASVAB prior to induction into the Navy. The 
test measures arithmetic reasoning, word knowledge, general 
science, mechanical comprehension, and other skills. On the 
basis of the ASVAB scores, recruits are divided into eight 


mental groups, based on percentiles as indicated in Figure 1. 


RECRUIT QUALITY MATRIX 


Parcentile Mental HIGH SCHOOL NON HIGH SCHOOL 
Group Diploma Graduate Diptoma Graduate 
33 
$3 


Armed Forces 
Qualification 
Test 


Figure 1: Recruit Quality Classification Matrix 








Figure 1 displays the complete classification matrix of 


prospective recruits. The Navy regards a prospect as a 
quality recruit if the individual scores in the upper 5Cth 
percentile on the ASVAB test, categories I to IIIU, and has or 
is expected to have a high school diploma. Quality recruits 


are also referred to as A-Cell contracts. 


B. FORECASTING A-CELL CONTRACTS 

In a report by Bohn and Schmitz [Ref. 1], NRC 
developed four statistical models to forecast the number of A- 
Cell contracts in each zip code area of CONUS. These models 


include the following variables as predictors: 


1. The target population of 17-21 year old males in the zip 
code: NRC estimated the population from the Navy’s 
Standardized Territory Evaluation and Analysis for 
Management (STEAM) database for 1991. 


2. The recruiter share assigned to each zip code: NRC 
assumes that a fractional number of recruiters or a 
recruiter share is assigned to each zip code. This 


fraction represents the amount of time recruiters devote 
to a particular zip code. 

The optimization models in this thesis allocate 
recruiter shares to zip codes in order to maximize total 
production of A-Cells. In the statistical models below, 
NRC simply assumes that recruiter shares are proportional 
to the 1991 population of 17-21 year old males since no 
record of recruiter shares was maintained for each zip 
code. 


3. The distance from the centroid of the zip code to its 
affiliated recruiting station: NRC used the longitude 
and latitude of the zip code centroids from commercially 
available software called MAPINFO [Ref. 2] to 
calculate the distances. In addition, NRC also assumes 
that the longitude and latitude of a recruiting station 





are those of the centroid of the zip code in which the 
station is located. 


4. The 1991 population density of the zip code: These 
densities are calculated from the physical zip code areas 
available from Litton Computer Services. 

The following describes the four forecasting models 
developed by NRC. The first three are log-linear and the last 


is quadratic. 


Log-Linear Regression Models 


Model 1: 
Acells, = aP'R§T!p$ = with R? = 0.53033 
where 
P, = the population of 17-21 year old males in zip code z 
R, = the recruiter share assigned to zip code z 
T, = the distance in miles from zip code z to its 
affiliated station 
D, = the population density in zip code z 


a = 0.08727, w = 0.64179, B = 0.28405, y = -0.32231, 


6 = -0.11923. 
Model 2: 
Acells,=aPzRitip? JJ e**! with R? = 0.53229 
1*3,5,7,8 
where 


Ss; = re if station is in Area pe 3c by a, B 


0, otherwise 
a@ = 0.08401, @ = 0.64506, 8B = 0.27904, y = -0.31380, 
6 = -0.11306, 





ky= 0.83109, ks= 1.03919, k,»= 0.96295, k= 1.15565. 


Model 3: 


6 
Acells, = aP¢RiT!p? J] e“!*! with R? = 0.53063 


i=2 


if #recruiters =i 


otherwise no 


if #recruiters 2 6 

otherwise 
a 0.05637, a@ 0.69148, B = 0.23419, y = -0.33103, 
6 = -0.12732 
k,=1.05817, k,= 1.06999, k= 1.15118, k= 1.09317, 


k= 1.17982. 


ic R ion Model 
Model 4: 


Acells, = a+ ®R,+yT,+8D,+ eT; +&R,T, with R? = 0.60598 


where 
a 0.51845, aw = 0.00160, B = 3.51236, y = -0.01095, 


6 -0.00172, € = 4.66E-5, £ = -0.01798. 


CONCAVITY ANALYSIS OF NRC PRODUCTION FUNCTIONS 

In economic theory, mathematical functions that relate the 
amount of inputs into a process to a specific level of output 
are known as econometric production functions. Clearly the 


NRC models are such functions, since they relate the expected 





number of A-Cell contracts produced to given amounts of 
inputs, such as recruiter share and target market population. 
An important class of production functions is the Cobb-Douglas 


function (Ref. 3] which has the form 


EUV x as ik a (2.1) 


where the function value f(v,,--,v,) represents the output of 
the process due to the nm inputs denoted by v,,--,v,. By 
definition, the v,’s are generally assumed to be nonnegative. 
However, when a, is negative, v, must be positive. To ensu-e 
that nonnegative inputs yield nonnegative outputs, A is 
generally a positive constant. 


The Cobb-Douglas function is homogeneous of degree 


k = @)+@,+-+@, [Ref. 4] because 


ECEV is SOV eo EV) = AEC Vee a) NB 02 1222) 


n 
In classical applications, the Cobb-Douglas functions also 
Satisfies the conditions that (i) k < 1 and (ii) a@ 20 for 
all ai = 1, ., om (Ref. 5]. This implies that the 
function ig decreasing return to scale, i.e., doubling all 
inputs results in less than doubled output. Furthermore, it 
can be shown that a Cobb-Douglas function satisfying the two 
conditions is strictly concave over the region where v, 2 0, 
i i, «;, 2 {[Reft: 6). The concavity plays an 
important role in optimization problems, for it guarantees 


global optimality. 








The first regression model considered by NRC, Model 1, is 


Similar in form to the Cobb-Douglas function. By redefining 
some variables, the function in Model 1 can be transformed 
into a Cobb-Dovglas function in terms of the variables R 
(recruiter share) and T (distance from each zip code to its 
affiliated station). In particular, the function is rewritten 


as 


£(R, T) = cRb Tr (2.3) 


where c is the positive constant aPD’. As stated, f(R,T) is 
not concave in R and T because y is negative. However, by 


defining w as 1/T, f(R,T) can be written as 


F(R, w) = cRPw (2.4) 


where } = -y. Since B+ A < 1 for Model 1, f(R,w) is strictly 
concave. Similar analysis can be performed for Models 2 and 
3 as well. 

In the above transformation, f(R,w) is undefined when 
T = 0. To avoid the problem of zero division, NRC assumes 
that the minimum distance between a zip code and its 
affiliated station is 0.5 miles, i.e., T 20.5. 


Model 4 is based on a quadratic function whose Hessian 


with respect to R and fT is given by 


H (x) -( ah (2.5) 


€ 2€ 





The eigenvalues of H(x) are 


0, = VE ve +e and 6,=-\F rere (2.6) 


Since ¢ = 4.66E-5 and — = -0.01798, @, > 0 and G5°<) Oy 2ee; 


H(x) is an indefinite Matrix. 


Therefore, Model 4 does not 


yield a concave production function. [Ref. 7] 
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III. OPTIMAL STATION ALIGNMENT FOR ONE RECRUITING DISTRICT 


The problem of aligning Navy recruiting stations within a 
single district involves selecting (or locating) recruiting 
stations to remain open and allocating recruiters to each of 
the open stations. The objective is to maximize the A-Cell 
contract prcduction. A preliminary discussion of this problem 
was provided in an interim report to NRC by Lawphongpanich, 


Rosenthal and Schwartz [Ref. 8]. 


A. RELATED RESEARCH 

Doll (Ref. 9] recently solved the problem. of 
locating Marine Corps recruiting stations. His objective is 
to maximize potential Marine Corps accessions as measured by 
a factor called the Propensity Weighted Qualified Military 
Available. In addition, Doll addressed the problem at the 
county level; i.e., he assumed that CONUS is a collection of 
counties, and the decision is whether or not to locate a 
recruiting station in a given county. Ina somewhat different 
problem, Celski [Ref. 10] used an uncapacitated 
facility location model to determine the optimal location of 
Army recruiting company headquarters. Celski'’s objective was 
to centrally locate a specified number of headquarters with 


respect to distance and population density of the target 


13 





market while simultaneously maximizing equity of recruiting 


responsibility. 


B. LOCATION-ALLOCATION PROBLEM FORMULATION 

To solve the station alignment problem in a given 
district, there are four sets of decisions to be made. The 
first set of decisions is to choose a subset of recruiting 
stations to open from a list of candidate stations. These 
candidates include existing and proposed stations. The second 
is to select the territory of each open station by assigning 
zip codes to stations. The third is to determine the number 
of recruiters to assign to each of the open stations. 
Finally, the fourth is to decide how to distribute the 
recruiter effort at a given station to zip codes in its 
territory. The main constraints are the number of stations to 
remain open and the number of available recruiters. The 
objective in aligning the stations is to maximize total A-Cell 
contract production in th. district which can be predicted by 
the first log-linear regression model in Chapter II. (The 
formulations involving the second and third log-linear models 
are Similar.) The mathematical formulation of the station 
alignment problem is called the location-allocation (LOCAL) 


problem and it is stated below. 


14 





C. MATHEMATICAL FORMULATION AND DISCUSSION 


INDICES: 
a = candidate station 
Zz = zip code under consideration 
DATA: 
T, = distance from zip code z to station i 
P, = population of 17-21 year old males in zip code z 
D, = population density of zip code z 
NR = number of recruiters in the district 
NS = number of stations to remain open in the district 
VARIABLES : 
X, = binary variable (=1 if station i is open and 0 
otherwise) 
Y, = binary variable (=1 if zip code z is assigned to 
station i and 0 otherwise) 
R, = number of recruiters to be assigned to station i 
SH, = proportion of recruiters (or recruiter. share) 


assigned to zip code z 








THE LOCATION-ALLOCATION (LOCAL) PROBLEM 
OBJECTIVE: 


y,,\* 
MAXIMIZE a+ P* » sy8 +( = * p> 
SH,R,X_Y, » af : » ioe ; 
CONSTRAINTS : 
YS xX, = Ns (Bek) 
1 
Vos bs da ea (332) 
> Veg sil, Vi (3.3) 
i 
> (Y,, * SH,) s Rj, Vi (3.4) 
Zz 
YR, < NR (3.5) 
i 


In the above formulation, parameters a, a, B, 6 and A are 
as defined in Chapter II. The objective is to maximize the 
total A-Cell contract production. Constraint (3.1) insures 
that the number of open stations equals the desired number, 
represented by NS. Constraint (3.2) requires zip codes to be 
assigned to open stations. Constraint (3.3) allows a zip code 
to belong to at most one station. However, the inequality 
allows unproductive zip codes to be left unassigned. 
Constraint (3.4) requires the recruiter shares distributed 
among the zip codes in the territory of station iI to not 
exceed the number of recruiters allocated to that station. 


Finally, constraint (3.5) guarantees that the number of 
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recruiters assigned to open stations does not surpass the 


number of recruiters, NR, available for the district. 

Since the objective function and constraint (3.4) are 
nonlinear, LOCAL is a nonlinear mixed integer programming 
problem which is not a well-solved problem in operations 
research [Ref. 11]. To solve LOCAL optimally would 
require a special algorithm to handle both the conditions of 
nonlinearity and integrality. However, the next chapter 
describes an alternative solution procedure which produces 
near optimal solutions efficiently and utilizes currently 


available software. 
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IV. A DECOMPOSITION APPROACH FOR THE LOCATION~ALLOCATION 


PROBLEM 


To avoid addressing the integrality and nonlinearity 
Simultaneously in LOCAL, the approach taken in this chapter 
decomposes the problem into four subproblems, which are solved 
in sequence. The solution to one subproblem is used as input 
to the subsequent subproblem. The following are brief 
descriptions of the subproblems. The sections below describes 


each of them in detail. 


1. The station location  subproblem: This subproblem 
determines the stations to remain open. 


2. The recruiter share allocation  subproblem: This 
subproblem assigns (possibly non-integer) numbers of 
recruiters to the open stations as determined by the 
first subproblem. 


3. The  integerization subproblem: When the second 
subproblem produces non-integer allocations of recruiters 
to stations, this subproblem optimally rounds those 
allocations to integer values. 


4. The recruiter share reallocation subproblem: Using the 
allocated recruiters provided by the third subproblem, 
this subproblem optimally proportions recruiter effort 
(or recruiter share) to zip codes in the territory of 
each open station. 

A. THE STATION-LOCATION SUBPROBLEM 
This subproblem assumes that recruiter shares have already 


been allocated to each zip code, i.e., the variable SH, in 
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LOCAL is fixed to some value, say SH, . Thus, the remaining 


decision variables are X, and Y,, which represent the decision 
to open or close station i and the assignment of zip code z to 
station i, respectively. This reduces LOCAL to the following 
station-location subproblen. 


STATION-LOCATION SUBPROBLEM 


OBJECTIVE: 
= y_.\4 
MAXIMIZE Yo a* Pp + SHE « » Z) * ps 
x,Y z 7 zi 
CONSTRAINTS: 
> xX, = Ns (4.1) 
i 
Voi 8 Aza Vizpa (4.2) 
> ¥og S 1, Vz (4.3) 
i 


As before, constraint (4.1) insures that the correct number of 
stations are opened, constraint (4.2) permits a zip code to be 
assigned only to open stations, and constraint (4.3) 


guarantees that a zip code is assigned to at most one station. 


B. THE RECRUITER SHARE ALLOCATION SUBPROBLEM 

An optimal solution to the station location subproblem, 
(X,Y), specifies the stations to remain open and the 
assignment of zip codes to stations. This provides a new 


alignment of recruiting stations within the district based 
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upon previously fixed levels of recruiter shares, SH, 


However, this allocation of recruiter share may not be optimal 
or feasible under the new alignment. For example, it would 
be infeasible for a zip code that was not assigned to any 
station to be allocated a recruiter share. The following 
nonlinear program optimally allocates recruiter shares to zip 
codes in the territory of an open station using the station 


alignment given by X and Y’, 


RECRUITER SHARE ALLOCATION SUBPROBLEM 


OBJECTIVE: 
MAXIMIZE )\ a * Pf + SHE « wt « Ds 
SH z 
CONSTRAINTS: 
DD Ya * SH, s NR (4.4) 
di Zz 


In the objective function, w, = Y, /T, where T, is defined to 
be distance from zip code zto station i. Constraint (4.4) 
prevents the recruiter shares distributed among zip codes from 
exceeding the available recruiters. The subproblem is an easy 
nonlinear programming problem since the recruiter share 
variable, SH,, is allowed to vary continuously for all z and 
there is only one constraint. The abjective function is 
strictly concave because 8 < 1. This guarantees solutions to 


the subproblem to be globally optimal. 
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C. THE INTEGERIZATION SUBPROBLEM 
Let SH, denote an optimal solution to the recruiter share 
allocation subp:oblem in the previous section. Then, the 


number of recruiters to be assigned to station i is: 


R,*> > ¥2; SH? Wi. (4.5) 


Zz 


These R, values may not be integer. Simple rules for rounding 
are insufficient for they may not be feasible. One feasible 


method of rounding R, is due to Rosenthal ([Ref. 12]. 


He first assigns R, recruiters to station i, where 


R, = FLOOR(R;). (4.6) 


This means that R,,, defined below, is the number 
Rees = NR -~ DUR, (4.7) 
i 


of recruiters left unassigned. Then, Rosenthal’s method 
allocates these residual recruiters to stations so as to 


maximize their marginal contract production. Note that R,, is 


integer since both NR and R, are integer. Since, R,, is 
expected to be much smaller than NR for practical problems, 


the degradation of the objective function from the continuous 


solution is very small. 
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The marginal productivity from adding u, recruiters to 


Station i is the difference between the expected A-Cell 


contract production with (R; + u,) recruiters and the 


production with only R, recruiters. The form of the NRC 


regression function (see Model 1 on page 8) does not permit a 
direct calculation of the marginal productivity, so the 
following approximation is used instead. 

Define WV, to be the set of all zip codes assigned to 
station i. Then, using the same notation as before, the 
productivity of station i is: 

0,= > ¢, su? 


oeee (4.8) 
where c, = aPSw*p. 


Now suppose the recruiter share, SH,, allocated to zip code z 
is proportional to the ratio of the target population in each 
zip code over the entire population in its affiliated 
station’s territory. Then, the recruiter share for zip code 


Z is defined to be: 
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P 
SH, = R; =". - forallzey,;. 


* SS P,, 


z'€y, 


(4.9) 


Using equation (4.9), the productivity for station i can be 


rewritten as 


0:= Ve (r ans ee , 
OE Be (4.10) 
z/éy, 


or equivalently 


) 
— pe _ P, 
Q,; = R; K where K, = c,|——*—_| . 
7 mse >» : 7") SS ep, (4.11) 


z/ ey, 


As defined, K, is constant for each zip code z. Thus, the 
marginal productivity for u, additional recruiters assigned to 
station i is: 


MO; = {(R, Pye Rf} > K, - (4.12) 


zey, 


Using this marginal productivity, the integerization 
subproblem can be formulated as follows: 


ADDITIONAL INDICES: 


3 = possible number of residual recruiters assigned toa 
station 
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ADDITIONAL VARIABLES: 


b, = binary variable (=1 if j residual recruiters are 
assigned to station i and 0 otherwise) 


INTEGERIZATION SUBPROBLEM 


OBJECTIVE: 
MAXIMIZE Yj} lb; * YK, * {(R, + 7)? - RP} 
b ig 264, 
CONSTRAINTS : 
Bee (F * Bis) S Rres (4.13) 
1 j 
YS bis, Vi (4.14) 


In constraint (4.13), the inner summation represents the 
number of additional recruiters to be assigned to station i 
and the constraint (4.13) itself insures that the additional 
recruiters assigned to all stations does not exceed the number 
of residual recruiters, R,,. Constraint, (4.14), guarantees 
that each station receives a unique number of additional 
recruiters, which can possibly be zero. Given an optimal 


solution, b,’, 
up = SoG * bi (4.15) 
j 


is the additional recruiters for station i, 


i.e., Rj = R; + uj. 





D. THE RECRUITER SHARE REALLOCATION SUBPROBLEM 

Solution to the integerization problem, R,, results in new 
recruiter allocations to each station i when the previous 
allocations were noninteger. Thus, the distribution of 
recruiter shares from the recruiter share allocation 
subproblem may not match the number of recruiters at each 
Station after the integerization. So, the recruiter share 
reallocation subproblem is to optimally reallocate recruiter 
shares to each covered zip code, given an integer assignment 
of recruiters to stations. Mathematically, this problem is 
formulated with the same variables and objective function as 
the recruiter share allocation subproblem, but constraint 


(4.4) is replaced by the following: 


YS (¥2, * SH,) s Ri Vi. (4.16) 


Zz 


This prevents the recruiter shares distributed to zip codes in 
the territory of station i from exceeding the number of 
recruiters allocated to that station as specified by the 


integerization subproblem. 





Vv. RESULTS FROM THE DECOMPOSITION APPROACH 


The principal focus of this chapter is to present the 
results for the Boston Recruiting District using the 
decomposition approach. However, first it is necessary to 
develop a bounding technique which can verify that the 


approach produces solutions with acceptable quality. 


A. BOUNDING THE DECOMPOSITION APPROACH 

Solutions obtained from the decomposition process are 
guaranteed to be both integer and feasible. To demonstrate 
the quality of the decomposition approach, its solutions are 
compared against a known upper bound to the LOCAL problem. 
Such a bound is determined by optimally solving a nonlinear 
programming relaxation (NLPR) of the LOCAL problem. 

The NLPR is obtained by omitting all binary requirements, 
i.e., the binary variables are allowed to range continuously 
in the interval [0,1]. Mathematically, the problem is 
formulated in the same manner as the LOCAL problem with the 
exceptions that constraints (3.4) and (3.5) are replaced by 
the following linear constraints: 


> SH, < NR (5.1) 


2 


SH, <3 7 *> Voie Vz. (5.2) 


i 
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Constraint (5.1) insures the recruiter shares distributed to 


zip codes do not exceed the available number of recruiters. 
Constraint (5.2) guarantees recruiter shares are distributed 
only to those zip codes assigned to stations; the parameter n 
is chosen to be large enough so that sufficient recruiter 
share can be assigned to a given zip code. With no binary 
restrictions, all of the decision variables are continuous 
and, as shown in Chapter II, the objective function is 
strictly concave. Consequently, the solution to the NLPR is 
unique and globally optimal [Ref. 13]. 

Table 1, on the following page, compares the results from 
the decomposition process against the corresponding upper 
bound for twenty-four problems using the 1991 recruiting data. 
Each row of the table specifies the district, number of 
stations and recruiters for the problem. Additionally, the 
values of recruiters and stations as a percentage of the 1991 
alignment is given within the parentheses. The percentage by 
which solutions from the decomposition approach differ from 
the upper bound is displayed in the last column. 

Table 1 shows that the decomposition erocess 
underestimates the true optimal solutions by at most 5.19 
percent and on average 2.5 percent. This difference is 
acceptable since the A-Cell production is predicted by a 
regression model with R’= 0.5303. Using the Amdahl 5990-500 


Computer at Naval Postgraduate School, the average solution 





TABLE 1: QUALITY OF SOLUTIONS FROM THE DECOMPOSITION 
APPROACH 









Total A-Cells 


[recruiters | _wurr | vecomp. | ¥ viteerence | 
jé_ta3.ay) [a3 (24a | 253.6: [ 2azas| 28 | 
101 














| aoe | 9 (23.7%) | 29 (25.0% | 420.60 
| 10a | 29 (50.08) | 58 (50.0%) | 570.78 


| 202 | 20 (74.28) _| 54 (74.08) 464.59 


27 (1008) | 73 (100%) | __ 535.54 
28 (73.7%) | 87 (75.0%) 674.94 | 672.09 fgg | 
38 (100%) | 116 (100%) | 758.53 


14 (50.08) | 52 (50.08) 
2a_ (75.08) | 78 (75.08) | 495.57] oes | 





2a (100%) | 104 (2008) | sea.62| sea.ze| 0.05 | 


7 (25.0%) | 26 (25.0%) 297.52 | 282.09 
6 (23.18) | 20 (25.08) 233.85 


161 

161 

161 

161 

101 

101 

101 

104 

104 

104 

104 

103 

103 

103 

103 

106 . . 
| 206 | 26 (200% | 80 roosy | ase.65| 496.63] 0.00 | 
| ao [6 (22.28) [29 aes) | 280.86] 273.32] 2.68 


27_(200%)_| 79_ (2008) 





28 


time for the twenty-four problems using the decomposition 
process is approximately 286 CPU seconds, while for the NLPR 
the time is approximately 19,275 CPU seconds. Thus, not only 
does the decomposition process produce near optimal integer 
solutions, but the solutions are obtained reasonably fast. 
Further analysis of Table 1 reveals that the numbers of 
recruiters and stations for which the decomposition process 
produces acceptable solutions ranges from approximately 25 to 
100 percent of the 1991 alignment. Moreover, the solution 
quality of the decomposition approach improves as the numbers 
of recruiters and stations are near 100 percent of the 1991 
alignment. To demonstrate, when the number of recruiters and 
stations are approximately 25 percent of the 1991 level, the 
average difference is 3.45 percent. When they are 
approximately 50, 75 and 100 percent of the 1991 level, the 
average differences are 1.92 percent, 0.73 percent and 0.12 
percent, respectively. This is due in part to the fact that 
the relaxed model, NLPR, yields more variables with integer 
values for problems with larger numbers of recruiters and 


Stations. 


B. RESULTS AND ANALYSIS 


The four optimization problems were implemented in the 


General Algebraic Modeling System [Ref. 14] (GAMS ) 


and the resulting integer and nonlinear programming problems 


were solved by the xS (Ref. 15] and MINOS 





{Ref. 16] solvers, respectively. The 1991 recruiting 
data for the Boston district is used to demonstrate 
applications of the model. For this example, the list of 
candidate stations consist of only those that were operational 
in 1991. 

In 1991, the Boston district had 42 recruiting stations 
and 111 recruiters. Figure 2 depicts results from the 
decomposition process with the number of stations varied from 
10 to 42 and the number of recruiters from 27 to 111. The 
graphs in this figure show how the number of A-Cells increase 


as a function of the number of stations and recruiters. 


Predicted Number of A-Cell Contracts for 
Boston Recrulting District 
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Figure 2: Predicted Number of A-Cell Contracts for Various 
Numbers of Recruiters and Stations for Boston District 


However, the expected gain due to the increase in the number 
of stations is decreasing. For example, the gain in the 
number of A-Cells by increasing the number of stations from 16 


to 25 is greater than that from either 25 to 33 or 33 to 42. 
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Information in Figure 2 can be used to construct isoquant 
curves representing combinations of stations and recruiters 
that yield the same number of A-Cell contracts. For example, 
the points where the horizontal line in Figure 2 intersects 
the five graphs represent combinations of recruiters and 
stations that produce an expected 600 A-Cell contracts for the 
Boston district. When several isoquants are constructed, they 
can be graphically displayed as in Figure 3. 


lsoquants for Boston District as a Percentage of the 1991 
Predicted Number of A-Cail Contracts 


Data cs 


757.01 (100%) 


881.31 (90%) 








Figure 3: Isoquants of Predicted Number of A-Cell 
Contracts for Boston District and Isocost Line 


To illustrate the use of the isoquants, suppose that the 
annual operating budget for Boston district, the average cost 
per recruiter and the average cost por staticn are $1,000K, 
$7.93K and 21.47K, respectively. Then the isocost line 
specified by the equation, 

7.93 x (mo. of recruiters) + 21.47 x (no. of stations) = 1000, 


represents all possible combinations of stations and 
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recruiters that satisfy the $1,000K budget exactly. This 


isocost line corresponds to the straight line shcwn in Figure 
a Since this iine is tangential to the isoquant for 
approximately 605.61 A-Cells, this is the maximum expected 
number of A-Cells obtainable with $1,000K. To obtain this 
number of A-Cells, Figure 3 shows that 20 stations and 72 
recruiters are necessary. On the following page, Table 2 
displays the stations that should remain open with 72 
recruiters. Each row contains the station identification 
number and the number of recruiters that should be assigned to 
that station. Note the total A-Ce.ls in Table 2 is 609.03 
instead of 605.61. This difference is due to. the 
interpolation used in constructing the isoquants in Figure 3. 

By varying the annual operating budget, different isocost 
lines can be drawn. Each would be tangent to different 
isoquants thereby yielding different combinations of 
recruiters and stations for the different budget limits. [In 
this manner, one can easily obtain the optimal numbers of 


recruiters and stations for any operating budget. 
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TABLE 2: RECOMMENDED ALLOCATION OF 72 RECRUITERS IN BOSTON 
DISTRICT WITH 20 STATIONS REMAINING OPEN 


Station I.D. No. Recruiters Expected 
a-Cells 


J 33.35 









31.00 


41.71 


34.11 


3 23.65 
28.02 
a ee 32.30 


609.03 


ie SO na 

ae 

ier eee Ee en Me 
a er eee ae ee 
eT ee ee Se 
ee er ee er eee 
er eee eee ee, 
peo Pca | 
ae aa a eee ee Ts 
se ee 


VI. ALLOCATING RECRUITERS AND STATIONS AMONG DISTRICTS 


The problem addressed in this chapter assumes that the 
current operating budget can support given numbers of stations 
and recruiters. These stations and recruiters are to be 
distributed among a collection of districts, such as a 
recruiting area. Intuitively, districts that produce more A- 
Cell contracts should receive more stations and recruiters. 
The overall objective is to maximize the total A-Cell contract 
production for the collection of districts. This problem is 
referred to as the multidistrict allocation (MDA) problem and 
can be approached in two ways. One possibility is to consider 
the collection of two or more districts as one ‘super’ 
district and use the method described in previous chapters to 
solve it. The other is to formulate the problem as a 
nonlinear resource allocation problem. 

The first approach has the disadvantage that the number of 
variables would be too large for solution in a reasonable 
time. For example, Recruiting Area 1, consisting of the 
following districts: Albany, Boston, Buffalo, New York, 
Harrisburg, Philadelphia, Pittsburgh and New Jersey, contained 
248 stations and 5,458 zip codes in 1991. Countirg the binary 
decision variables alone, the problem of allocating recruiters 
and stations to districts in Area 1 can contain up to 248 + 


(248 x 5,458) or 1,353,832 variables. Thus, the resource 
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allocation formulation, which has a much smaller number of 


variables, is chosen instead. The difficulty with this 


formulation is that the objective function is not expressible 
in closed form. In one of the sections below, the objective 
function is approximated using piecewise linear functions. 
Oth.sr sections formally state the formulation, describe an 


implementation and present an example. 


A. THE MULTIDISTRICT RESOURCE ALLOCATION FORMULATION 

In this allocation problem, a number of recruiters and 
stations are to be distributed to two or more districts so as 
to maximize the total A-Cell production. Given that there are 
costs associated with employing a recruiter and maintaining a 
station, there is also a constraint to limit the cost not to 
exceed a certain budget. Below is a formulation of this 
problem. 


MULTIDISTRICT ALLOCATION (MDA) PROBLEM 
INDICES: 


d district 


total annual budget for the districts considered 
number of recruiters to be allocated 
number of stations to be allocated 


average annual cost for a recruiter ir district d 








cS; = average annual cost of a station in district d 


DECISION VARIABLES: 


Ke number of recruiters assigned to district d 


Sa number of stations assigned to district d 
OBJECTIVE: 


MAXIMIZE > Fy (Rg, Sq) 
R,S d 


CONSTRAINTS : 


> Ry < NR (6.1) 
d 
¥ Sq < NS (6.2) 
d 

Yo (ciRy + CGSq) s B (6.3) 


d 

In the objective, F,(r,s) represents the maximum A-Cell 
production using r recruiters and s stations in district d. 
Since F,(r,s) is generally nonlinear, the above problem is a 
nonlinear integer programming problem. Its first two 
constraints, (6.1) and (6.2), insure that the number of 
assigned recruiters and stations do not exceed the available 
number and the last, constraint (6.3), ensures that the cost 


does not exceed the budget. 
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The function F,(r,s) is not expressible in closed form. 
To evaluate F,(r,s) at a particular value of r and gs, the 
LOCAL problem corresponding to district d must be solved. As 
seen in Chapter V, when all the variables in LOCAL are allowed 
to vary continuously, the resulting problem is a strictly 
concave programming problem. The function F,(r,s) is known as 
an optimal value function and can be shown to be concave (Ref. 
13]. When some of the variables in LOCAL are restricted to 
integers, the concavity for F,(r,s) is not guaranteed. 
However, our experiments empirically support the assumption 
that F,(r,s) is concave when s is integer and r is allowed to 
range continuously. 

Based on the magnitudes of R, and S,, the integer 
restriction of R, may be relaxed, but not S,. Using this form 
of relaxation, F,(r,s) can be approximated with piecewise 
linear functions and the MDA problem can be reformulated as a 
linear integer programming problem. The next section 


describes the approximation scheme and the linear formulation. 


B. APPROXIMATING THE MULTIDISTRICT ALLOCATION PROBLEM 

Let f£,,(r} denote the value of F,(r,s) at a fixed value of 
s. Since it is assumed that noninteger recruiters are 
allowed, f,(r) is a continuous function of r and can be 
approximated with a piecewise linear function g,,(r). 

By virtue of the assumption that +,(r) is piecewise linear 


and concave, its value can be maximized by linear programming 
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{Ref. 17]. As an illustration, suppose g,,(r) is 
specified by the points {u,,g,,(u,)} for i=1,...,4 (see Figure 
4) where g,(u,) represents an actual solution to the LOCAL 
problem with u, recruiters allocated to s stations in district 


d. 


22) 
B 
a 
s 
[sj 
o 
is) 
2 
° 
< 
2 
o 
5 
3 
2 
a 


Ue 
Recruiters 





Figure 4: A Piecewise Linear Concave Function 


Let p, be the proportion of the recruiters r represented in the 
interval between u,,; and u, where u, > U,,;. Then, any 


Os rs uy can be written as 


re py Vp, Sup up, (6.4) 
iet 
where p, i8 nonnegative and p,,, > 0 only when p, = u, For 


example, let {u,, u,, Us, Us} = {10, 20, 30, 40}, then for 
r= 25, {P,, Do, Dse De} = {10, 10, 5, 0}. Thus, by defining r 


in this manner, g,,(r) can be written as 
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4 
9a,s(2) => mp, {6.5) 


ix1l 

where m represents the slope of the piecewise linear segment 
between u,, and u,. 

In the formulation below, the approximating function g,, (r) 
is not computed for every possible value of s. Instead, g,,(r) 
is computed for only a small collection of s values. The 
remaining g,(r) are obtained by linear interpolation. For 
example, assume that g,,,(r) and g,,;;(r) have been computed, 


then g,,,(r) is approximated as follows: 


Ja,i2'F) = 394,10 (2) i 294,15 (7) « (6.6) 


Mathematically, the approximation to the MDA problem is 


formulated below. 


APPROXIMATE MULTIDISTRICT ALLOCATION PROBLEM 


INDICES: 
d= district 
Ss = number of stations to remain open 
k = linear segment for the piecewise linear functions 
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GIVEN AND DERIVED DATA: 


cS, = average annual cost of a station in district d 
C average annual cost for a recruiter in district d 
B total annual budget 
M,, = the siope of segment k fcr the piecewise lineai 
function approximating f,, (r) 
U, = upper bound for the interval of segment k for 
district d 
SU, upper limit on the number of stations which may 
remain open within district d 
Sly lower limit on the number of stations which may 
remain open within district d 
NS number of stations available 
NR number of recruiters available 
VARIABLES : 
Pak =  mnumber of recruiters in the interval for segment k 
with s stations open ia district ad 
q, = binary variable(=1 if s stations are to remain open 
in district d and 0 otherwise) 
Va slack variable 
OBJECTIVE: 
MAXIMIZE OSS Masc* Pace 
Pq dos k 
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CONSTRAT 


oo 


Pask S Yay * Vass Vd,s,k (6.6) 
s dag = 1, Vd (6.7) 

3 
yy (8 * aa) < NS (6.8) 

dad s 


pp dp > Pask S$ NR (6.9) 
s 


Vat (8 * dys) = SUg, Wad (6.10) 
3 
S*das5 > Pax Vas (6.11) 
k 


LLY (ease * Cr) + VY (das * Ca) s B (6.12) 
ds k dis r 


In the above formulation, constraint (6.6) ensures that the 
correct segments of the piecewise linear functions are used 
for each district. Constraint (6.7) ensures a unique number 
of stations remain open within each district. Constraint 
(6.8) allows at most NS stations to remain open. Constraint 
(6.9) insures the number of recruiters to be allocated, NR, do 
not exceed the number available. 

Constraint (6.10) is a two-sided constraint represented by 
a single constraint [Ref. 18]. It guarantees) that 


the number of stations open in each district is between a 
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lower and an upper bound. A lower bound is established to 
avoid potential political problems by not recruiting within 
specific districts or conversely recruiting exclusively in 
others. The upper bound prevents the number of stations 
within any district from exceeding the 1991 level, since NRC 
is downsizing. 

Constraint (6.11) guarantees that the number of recruiters 
assigned to any district is at least as large as the number of 
stations. Finally, constraint (6.12) insures the costs 
associated with the number of recruiters and stations 


allocated do not exceed the available annual budget. 


C. EXAMPLE 

To illustrate, the above formulation was implemented in 
GAMS and solved by OSL ([Ref. 19]. The data inputs 
are from Recruiting Area 1 which, as mentioned previously, 
consists of eight districts. In 1991, Area 1 had 248 station 
to which 702 recruiters were assigned. NRC estimated that the 
average annual costs for recruiters and stations for all 
districts in Area 1 are approximately $50K and $9.472K, 
respectively. Using these estimates, the 1991 operating 
budget for Area 1 is calculated as 

(248 x $9.472K) + (702 x $50K) =~ $37,449K. 

On the following page, Figure 5 depicts solutions to the IA 
problem for budget limits ranging from 50-90 percent of the 


1991 budget. To obtain the graph in this figure, a total of 
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five multidistrict problems must be solved which takes on the 
average 3.655 CPU seconds using the Amdahl Computer. The 


graph in this figure shows how the number of A-Cell contracts 


increase as a function of the operating budget. 


Estimated A-Cell Contracts 


220 240 260 280 300 
Budget Limit { Thousands of dollars) 


Figure 5: Solutions to the MDA problem for Area 1 


Based upon the graph in Figure 5, the increase in A-Cells 
appears to be fairly constant with respect to budget 
increases. This is because it is more cost effective to 
increase the number of recruiters instead of opening stations 
as the budget increases. 

In particular, consider a ten percent reduction in the 
budget, i.e., solve the multidistrict allocation problem for 
Area 1 with $33,699K instead of $37,449K. This yields a new 
alignment displayed in Table 4. When compared to the A-Cell 


contracts produced under the 1991 alignment (i.e., 4,432 





contracts), the new one yields more contracts by approximately 


two percent. 


TABLE 4: A SOLUTION TO THE MDA PROBLEM FOR AREA 1 


| Districts Number of Number of Estimated Cost Predicted 
Racruiters Stations (in thousands of $) A-Cell 
Contracts 





627 33699 .27 4522.19 
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VII. SUMMARY AND CONCLUSIONS 


This thesis addresses two problems in aligning the 
recruiting structure for Navy Recruiting Command. The first 
problem involves two decisions affecting recruiting stations 
within a single recruiting district: which stations should 
remain open and how many recruiters to assign to each open 
station? The second problem is to decide how many recruiters 
and stations each district should have. 

The first problem is modeled as a nonlinear mixed integer 
programming problem. To obtain a solution with readily 
available software, the problem is decomposed into four 
subproblems solved sequentially. They are (i) the station 
location, (11) the recruiter share allocation, (iii) the 
integerization and (iv) the recruiter share reallocation 
problems. This decomposition approach produces good integer 
solutious relatively fast. These solutions are within 2.5 
percent of optimality on the average and within 5.19 percent 
in the worst case. Results from the decomposition approach 
were used to construct isoquants for A-Cell production in the 
Boston district. When superimposed with an isocost line, 
these isoquants identify the most economical number of 
recruiters and stations for the Boston district. 

To keep the number of variables at a manageable size, the 


second problem is modeled as a nonlinear resource allocation 
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problem in which the objective function is not available in 


closed form. To efficiently solve this problem, the function 
is approximated in a piecewise linear fashion using the 
results from the first problem. The approximation yields a 
linear integer programming problem which can be solved by 
readily available software. For a given operating budget, 
this integer programming problem was used to produce optimal 
numbers of recruiters and stations to be allocated to each 
district in Recruiting Area 1. Using these results, it is 
shown that the 1991 A-Cell production level in Area 1 can be 
maintained using only 90 percent of the 1991 operating budget. 

In conclusion, this thesis provides NRC with a tool to aid 
in their current downsizing efforts. The effectiveness of 
this tool largely depends on two unique methodologies 
developed in earlier chapters. One utilizes decomposition to 
obtain solutions to nonlinear mixed integer programming 
problems for which no readily available algorithm exists. 
Unlike other heuristic approaches, the decomposition produces 
answers that can be considered optimal for all practical 
purposes. The other methodology uses an alternate formulation 
and an approximation scheme to solve problems that are 


otherwise too large for currently available solvers. 


A. AREAS FOR FURTHER RESEARCH 
The results of this thesis also point out the following 


areas for further research. 
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1. Robustness of the Model Against Different Forecasting 
Models for Contract Production 
This thesis uses the log-linear form of the 
forecasting model as recommended by NRC. However, there are 
other competing forecasting models which are of different 
forms and include other demographic and socio-economic 
factors. If different forecasting models generate drastically 
different station locations and recruiter allocations, then 
the choice of forecasting models will play an important role 
in the realignment decisions. However, if different 
forecasting models generate similar locations and allocations, 
then the forecasting model which expedites the solution 
process should be chosen. In any case, this discussion 
highlights the need for further investigation to determine 
model robustness. 
2. Realignment System Development 
It is inevitable that the need to realign the 
recruiting organization will re-occur in the future. This may 
be due to a change in policy or a shift in the location of the 
target market. To keep recruiting efficient, the recruiting 
organization should be responsive to these changes. This will 
require that the optimization models developed in this thesis 


be integrated in a user-friendly system. 
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